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Abstract. We review the constraints given by the linear matter power spectra data on cosmological and bias 
parameters, comparing the data from the PSCz survey (Hamilton et al., 2000) and from the matter power spectrum 
infered by the study of Lyman alpha spectra at z=2.72 (Croft et al. , 2000) . We consider flat -A cosmologies, allowing 
A, Ho and n to vary, and we also let the two ratio factors Vpscz and riyman {ri — pj^'j^(fc) ) vary independently. 
Using a simple minimisation technique, we find confidence intervals on our parameters for each dataset and 
for a combined analysis. Letting the 5 parameters vary freely gives almost no constraints on cosmology, but 
requirement of a universal ratio for both datasets implies unacceptably low values of Ho and A. Adding some 
reasonable priors on the cosmological parameters demonstrates that the power derived by the PSCz survey is 
higher by a factor ~ 1.75 compared to the power from the Lyman a forest survey. 
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1. Introduction 

Bias, defined to be the ratio of rms luminosity density 
fluctuations to dark matter fluctuations at a fiducial scale 
of 12/1^^3 Mpc, is the most uncertain of the minimal set of 
cosmological parameters. These may be taken to be A, Hq, 
rig and Qb for the background, with the inhomogeneities 
described by n and the bias parameter b. 

There are three techniques available for computing 
the shape and amplitude of the power spectrum of den- 
sity fluctuations. These involve using the CMB anisotropy 
data (BOOMERanG, MAXIMA, DASI), galaxy redshift 
surveys (PSCz, Hamilton & Tegmark, 2000) and the 
Lyman alpha forest (Croft et al., 2000). Use of large- 
scale peculiar velocity surveys and of galaxy cluster abun- 
dances provide a value of the normalisation, via the pa- 
rameter dg^^o'^, where erg = 6 evaluated at the fiducial 
scale of 12/i2/3 Mpc. We confirm that CMB and galaxy 
surveys provide a consistent data set for the power spec- 
trum. However incorporation of the matter power spectra 
inferred from Lyman alpha forest requires a relative bias 
of galaxies relative to the Lyman alpha forest once reason- 
able priors are adopted for the cosmological parameters. 

2. Data sets 
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2.1. PSCz 

The linear power spectrum inferred by the PSCz survey is 
taken from Hamilton & Tegmark 2000 (HTOO hereafter). 
We choose the 22 points given in table B2 of HTOO, which 
correspond to estimates of the decorrelated linear power 
spectrum. The dataset is shown in the upper part of figure 

B 

2.2. Lyman a 

We use the linear matter power spectrum inferred from the 
flux power spectrum of lyman a forest spectra by Croft et 
al., 2000 (COO hereafter). This corresponds to 13 estimates 
of the spectrum at < z >= 2.72, referenced in table 4 of 
COO. When these data are used, the normalisation uncer- 
tainty (-1-27%, -23%) is included, then marginalised over 
(see next paragraph). The data are drawn in the lower 
part of figure ^ 

3. Cosmologies 

3.1. Free parameters 

In this paper, we consider flat cosmologies, with no reion- 
isation, no massive neutrinos and no tensor fluctuations. 
We also fix the baryon fraction to be consistent with Big- 
Bang nucleosynthesis and CMB determinations: ilhh^ = 
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Fig. 1. Actual data set of matter power spectra. Up: PSCz 
spectrum. Down: Lyman alpha power spectrum. The solid 
line is the best model of the combined analysis (see text) 
and the (blue) dotted line correspond to the best model 
of the CMB, PSCz, Lyman a combined analysis. 

0.020. The remaining parameters are listed in table 1, as 
well as the corresponding ranges and steps. All the theo- 
retical power spectra are COBE-normalised, and the ratio 
parameters are defined as: 

P{k)pscz flyman — "^pscz /lyman^crnbik) ■ 

Here Ao is the vacuum density and n is the fluctuation 
spectral index. 





Ho kms~^/Mpc 


Ao 


n 




'^lyman 


Min. 


20 


0.0 


0.70 


0.50 


0.10 


Max. 


100 


1.0 


1.30 


1.475 


2.05 


step 


10 


0.1 


0.03 


0.025 


0.05 



Table 1. Parameter space explored: 

f^tot = 1 



3.2. Theoretical matter power spectrum 

We have considered theoretical matter power spectra with 
an arbitrary val ue of the sha pe p arameter (F), u sing 
the approach of Bardccn 1986 and Sugiyama 1995 . We 
computed t he COBE normalisation of the spe ctra using 
the CAMB ( [Lewis, Challinor, fc Lasenby 2000[) code. We 
used the "grow package" from Hamilton 2001 to derive 



the growth factor. Our two theoretical power spectra at 
z=0 and z— 2.7 are presented in the following equations: 



Pz=o{k) = NiX TF{V,kf X g{kf x Po{k) * r 



2 

pscz 

Pz=2.72{k) ^ N^xf^x TF{T,kf X g{Kf x P,{k) * rf^^^ 

where we write: 

= (l + z)-2*[(A + a„*(l. + z)3 + f],*(l + z)2)i/2]3(3) 

and the initial matter power spectrum is given by: 
Po(fc) - (fc * hT (4) 
We take the transfer function to be: 



(1) 



TF(r,fc).Mi±^:M^x 

^ ' ' (2.34 *g) 



(1 + a2 * 9 + ("3 * qY + ("4 * qY + (as * qYY'^ 



(5) 



where q = fc/F. We use the coefficients given by Sugiyama 
(1995): ai = 2.34, = 3.89, = 16.1, = 5.46, as = 
6.719. Here we treat F as a free parameter. We then 
marginalize over it. The shape of the power spectrum will 
be investigated in a future paper. iVi and N2 are normal- 
isation factors, such that our matter power spectra are 
COBE-normalised. 

4. Statistics 

4.1. Likelihood 

We approximate the likelihood functions of our param- 
eters by independent multivariate Gaussians, using the 
following expression: 



-'pscz flyman CXp 



^pscz / lynian 



For the PSCz data, the has been computed as follows 
(P,(fc^,)--P™°''''(0p|fc»))' 



Xpsczi^p) ^ ^ ] 



(6) 



where p">-°d.<^'- is derived from equation |l| of the following 
section and Op is the set of parameters explored listed in 
table 0. 

Concerning the Lyman alpha data, we take into con- 
sideration the uncertainty on the normalisation. We add 
a contribution to x^, assuming a double-tail Gaussian for 
the normalisation distribution: 



xlymani^liN) = X! 



Xlymani'^l) = -2 X log 



{p,{k,)-p{euN\hf 



yman 



)dN 



(7) 
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4.2. Results 



The results are shown in 2D contours plots, with red- 
dashed lines corresponding to confidence intervals 68, 95, 
99% in one dimension, and filled contours corresponding to 
68, 95, 99% in two dimensions. The remaining parameters 
have been marginalised (not integrated), finding the min- 
imum for each pair of parameters plotted in the graphs. 



5. Results and constraints 

5.1. PSCz and Lyman a matter spectra 

Given our 5 free parameters, the constraints given by each 
set independently are almost inexistent. Many combina- 
tions of parameters lead to degeneracies which allow values 
even out of our grid of parameters. 

5.2. Are the two data sets consistent 7 

We have seen that our two data sets allow most 
of the values for the parameters we are considering. 
Nevertheless, this does not mean that they are compat- 
ible. Marginalising over parameters, for a better view 
of the constraints, does not allow us to include all the 
correlations and degeneracies between parameters. This 
is even more critical when the constraints are weak. To 
see whether the data are consistent in some particular 
cosmologies, we adopt two scenario. The first is to con- 
sider that the ratios are the same for both datasets. 
Marginalising over these will reveal the prefered range for 
the cosmological parameters. The second approach is to 
let the ratios vary freely, while putting some priors on the 
cosmological parameters. This will give us estimates of the 
ratio for each data set. 

5.2.1. Uniform ratio 

We have combined the data sets by adding the two 
grids defined in equations ^ and ^ r being the unique 
bias parameter. The prefered models obtained with this 
assumption correspond to very low values of Ho and A. 
The "best one" is plotted in figure ||. The value of the 
goodness of fit for this model, given by the absolute value 
of the on the two dataset, is good {Xmin = 22 for 30 
degrees of freedorrj^ . The corresponding confidence inter- 
vals are shown in figure |^, where the remaining parameters 
have been marginalised over. This combined analysis give 
strong constraints on the cosmological parameters, which 
are largely in contradiction with most estimates of the cos- 
mological parameters from independent techniques (such 
as SNIa, age of the Unviverse, CMB). Although the con- 



^ The degrees of freedom are given by the number of points 
in the minus the number of free parameters investigated; 



dof = Npts - Np. 



This is actually correct only if the x is 
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Fig. 2. Combined constraints from the Lyman a power 
spectrum and the PSCz spectrum with the same ratio. 
Filled contours mark the 68, 95 and 99% confidence levels 
in 2-dimensional space. The dashed red line define the 68, 
95 and 99% confidence levels when projected on 1 param- 
eter. 

cordance model is marginally acceptable, a low tilted 
CDM model is preferred. 

5.2.2. Ratio parameters 

The preceding discussion does not reinforce our faith in 
the concordance model with respect to other parameter 
estimates. The analysis of the data could be done via an al- 
ternative approach. Assuming that we "know" the cosmo- 
logical parameters, we may examine the relative bias of the 
two data sets. Statistically speaking, this means that we 
can put priors on the cosmological parameters, marginalise 
over them, and see what the constraints are on the 2 re- 
maining parameters that are associated with the inhomo- 
geneities power. We adopt the reasonable following priors 
in our analysis: A = 0.7±0.1, iJo = 65±10 kms~^Mpc~^ , 
n = 1 ± 0.1 which are in agreement with most of the ac- 
tual parameter estimates in the literature. We assume a 
Gaussian shape for the priors, giving the following likeli- 
hood: 



linearly dependent on the parameters, but it is generally used, 
and still gives a good estimate of the goodness of fit. 



X'(A) 
X^(Ho) 
X^in) 



Xlscz + Xfyman + X^W + X^(^o) + x"^ (n) 

(A -0.7)' 
0.12 

{Ho - 65)2 

102 
(n- 1.0)2 

(112 
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Fig. 3. Combined Constraints from the Lyman a matter 
power spectrum and the PSCz spectrum with priors on 
cosmological parameters (see text) 



Fig. 4. CMB constraints using COBE, BOOMERanG, 
DASI and MAXIMA datasets for flat ^^h^ = 0.02 models. 



then 



")(.corahirLe 



di'>'pscz,riy,nan) = Min[A.H„.n]iX ) (8) 

Figure ^ shows the confidence intervals for Vpscz and 
flyman- The bcst modcl is obtained for Tp^cz = 1.075 and 
flyman = 0.6, with a good goodness of fit: = 24.8 for 29 
degrees of freedom. We can clearly see, in the figure I the 
correlation between the two biases that leads to a relation 
(over the explored range of parameters and between the 
95% confidence level interval): 



' pscz 



— 1.57 * n 



yman 



0.16 or ^^'"^ ~ 1.8 (9) 

flyman 

The corresponding 68% intervals on each ratio are: 0.45 < 
flyman < 0.8 and 0.8 < riynian < 1.35; the constraints 
being low at 95% CL. 

5.2.3. Combining CMB and matter spectra 

Instead of setting independent priors on each cosmologi- 
cal parameter, we can use the latest CMB analysis as a 
prior. In LeDour et al., 2000 and Douspis et al., 2001, we 
used the current CMB data set to derive constraints on 
cosmological parameters. The latter paper considers the 
most recent data including BOOMERanG (Netterfield et 
al, 2001), DASI (Halverson et al., 2001) and MAXIMA 
(Lee et al., 2001). We use the same grid as in table 1 
for the CMB analysis. As seen in the literature (Douspis 
2000, Jaffe et al., 2001, Melchiorri and Griffiths, 2000 and 
others), the CMB by itself does not constrain Hq or A; 
both of these parameters are degenerated with the total 
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Fig. 5. Combined Constraints from CMB, PSCz spectrum 
and the Lyman a power spectrum in the {r pscz, flyman) 
plane 



density (or the curvature). As we have fixed fifjft,^ = 0.02 
and ^totai = 1 in this analysis, the CMB now gives inter- 
esting confidence intervals on the remaining cosmological 
parameters (see figure ^). These constraints are similar 
to those used in the previous section, but contain more 
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information, because of the correlations between parame- 
ters. This way of setting priors should be preferable to the 
previous one. 

We constructed a new likelihood grid from the CMB 
and the matter power spectra: 



X {d) — A.cmb + Xpscz + Xlyman 

with B = {^^,H(i,n,T,rpscz,rlyynan) 



(10) 
(11) 



where A is the approximation developed in Bartlett et al., 
2000. We then maximise this new grid and marginalise to 
find confidence intervals. The best model given by A = 
0.6, Hq = 60, n = 0.97, Vpscz = 0.875, ny^an = 0.5 has 
a good goodness of fit and is plotted in figure |l|. 

As the constraints on cosmological parameters are 
stronger in the CMB analysis than in the matter power 
spectrum analysis, we do not expect any improvement on 
cosmological parameter confidence intervals. On the other 
hand, marginalising over A, Hq, n gives us interesting 
constraints on Tpscz and riyman- When the three datasets 
are combined, we find constraints on the relative bias and 
again a degenerated relation between the two ratio param- 
eters, leading to a lower ratio factor for the Lyman linear 
power spectrum (cf figure ^) . The 68 % confidence level are 
the following: 0.45 < riy^an < 0.6 and 0.7 < nyman < l-l- 
We may also express the relative bias between the two sets 
of data as: 



1.33 * riy„ 



0.20 or 



'^lyman 



1.75 



(12) 



lower bias, or even an antibias, compared to that of lu- 
minous galaxies. Dwarfs are likely to be the final fate of 
much of the forest. Indeed LCDM simulations of cluster- 
ing find that luminous galaxies themselves are somewhat 
antibiased. Hence our result fits in well with hierarchical 
structure formation, since the low Lyman alpha forest ra- 



tio 



1.75 is similar to that of dwarf galaxies. 



A., 



Le Dour, M. 
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6. Conclusions 

Determination of the cosmological model parameters is 
entering into a new era of precision. The weakest link in 
our understanding of the universe is in the realm of the 
density fiuctuations, and their evolution. It is well known 
that bias depends on the objects being sampled, and varies 
with galaxy luminosity and morphological type. However 
these variations are attributed to the complex physics of 
galaxy formation. We find here that by adopting reason- 
able cosmological model priors, the matter power spec- 
trum inferred by Lyman alpha forest power spectrum can 
be shown to be biased low relative to the the bias of 
galaxies. 

It is usually assumed that the Lyman alpha forest, rep- 
resenting gas that has not yet virialized and has density 
contrast between unity and 200, is a good tracer of the 
dark matter (e.g. Croft et al. 2001). In fact this is a dan- 
gerous assumption, since the Lyman alpha forest gas is a 
very small fraction of the total gas density, 99.99 percent 
or more of which is ionized. Moreover, the proximity effect 
observed for quasars and more recently for Lyman break 
galaxies (Steidel 2001), demonstrates that the forest is not 
a good tracer of the total matter density on scales of up 
to tens of megaparsecs. 

We note that gas-rich dwarf galaxies, which avoid the 
vicinities of luminous galaxies, also display a significantly 



